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Geometry and material effects in Casimir 
physics - Scattering theory 

Sahand Jamal Rahi, Thorsten Emig, and Robert L. Jaffe 



Abstract We give a comprehensive presentation of methods for calculating the 
Casimir force to arbitrary accuracy, for any number of objects, arbitrary shapes, 
susceptibility functions, and separations. The technique is applicable to objects im- 
mersed in media other than vacuum, to nonzero temperatures, and to spatial arrange- 
ments in which one object is enclosed in another. Our method combines each ob- 
ject's classical electromagnetic scattering amplitude with universal translation ma- 
trices, which convert between the bases used to calculate scattering for each object, 
but are otherwise independent of the details of the individual objects. This approach, 
which combines methods of statistical physics and scattering theory, is well suited to 
analyze many diverse phenomena. We illustrate its power and versatility by a num- 
ber of examples, which show how the interplay of geometry and material properties 
helps to understand and control Casimir forces. We also examine whether electro- 

t— I dynamic Casimir forces can lead to stable levitation. Neglecting permeabilities, we 

prove that any equilibrium position of objects subject to such forces is unstable if the 
permittivities of all objects are higher or lower than that of the enveloping medium; 

-y, the former being the generic case for ordinary materials in vacuum. 
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1 Introduction 

Neutral objects exert a force on one another through electromagnetic fields even 
if they do not possess permanent multipole moments. Materials that couple to 
the electromagnetic field alter the spectrum of the field's quantum and thermal 
fluctuations. The resulting change in energy depends on the relative positions 
of the objects, leading to a fluctuation-induced force, usually called the Casimir 
force. Alternatively, one can regard the cause of these forces to be spontaneous 
charges and currents, which fluctuate in and out of existence in the objects due 
to quantum mechanics. The name 'Van der Waals force' is sometimes used in- 
terchangeably but it usually refers to the Casimir force in the regime where ob- 
jects are close enough to one another that the speed of light is effectively infinite. 
The Casimir force has been the subject of precision experimental measurements 

E2|62|8a|23[l3|23|2i|4i[nD|52|2^ and can influence 

the operation of nanoscale devices |[T6l[T3l . see reference ll54l for a review of the 
experiments. 

Casimir and Polder calculated the fluctuation-induced force on a polarizable 
atom in front of a perfectly conducting plate and between two polarizable atoms, 
both to leading order at large separation, and obtained a simple result depending 
only on the atoms' static polarizabilities [ 15 ]. Casimir then extended this result to his 
famous calculation of the pressure on two perfectly conducting parallel plates fl4l . 
Feinberg and Sue her 0411351 generalized the result of Casimir and Polder to include 
both electric and magnetic polarizabilities. Lifshitz, Dzyaloshinskii, and Pitaevskii 
extended Casimir's result for parallel plates by incorporating nonzero temperature, 
permittivity, and permeability into a general formula for the pressure on two infinite 
half-spaces separated by a gap [64, 25l l65l . 

While these early theoretical predictions of the Casimir force applied only to 
infinite planar geometries (or atoms), the first precision experiments measured the 
force between a plate and a sphere. This geometry was preferred because keeping 
two plane surfaces parallel introduces additional challenges for the experimentalist. 
To compare the measurements with theory, however, a makeshift solution had to be 
used: known as the Proximity Force Approximation (PFA), it estimates the Casimir 
force by integrating the Casimir pressure between opposing infinitesmal surface area 
elements, as if they were parallel plates, over the area that the sphere and the plate 
expose to one another f74l . In general, this simple approximation does not capture 
curvature corrections but in many experimental situations, it performs surprisingly 
well, as can be seen in Fig.[T] for example; at the small separations at which the force 
is typically probed in precision measurements the sphere and the plate surfaces are 
well approximated by a collection of infinitesimal parallel plates. 

Clearly, for larger separations and for surfaces that are not smooth, the PFA must 
fail. For example, in measurements of the Casimir force between a sphere and a 
trench array significant discrepancies were found ifTTl . And even for the regimes 
in which the PFA yields good estimates it would be desirable to know what the 
corrections are. 
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Fig. 1 Force between a sphere of radius « 100/im and a plate, both coated with Au-Pd |69|. 
Square dots represent measurements, the solid line is a theoretical computation using the PFA 
approximation and taking into account roughness and finite temperature corrections as well as 
material properties. The other lines represent calculations, where some of these corrections are not 
taken into account. 



In order to study Casimir forces in more general geometries, it turns out to be 
advantageous to describe the influence of an arrangement of objects on the elec- 
tromagnetic field by the way they scatter electromagnetic waves. Here, we derive 
and apply a representation of the Casimir energy, first developed with various lim- 
itations in Refs. |28, 29] and then fully generalized in Ref. [75], that characterizes 
each object by its on-shell electromagnetic scattering amplitude. The separations 
and orientations of the objects are encoded in universal translation matrices, which 
describe how a solution to the source-free Maxwell's equations in the basis appro- 
priate to one object looks when expanded in the basis appropriate to another. The 
translation matrices depend on the displacement and orientation of coordinate sys- 
tems, but not on the nature of the objects themselves. The scattering amplitudes and 
translation matrices are then combined in a simple algorithm that allows efficient 
numerical and, in some cases, analytical calculations of Casimir forces and torques 
for a wide variety of geometries, materials, and external conditions. The formalism 
applies to a wide variety of circumstances, including: 

• n arbitrarily shaped objects, whose surfaces may be smooth or rough or may 
include edges and cusps; 

• objects with arbitrary linear electromagnetic response, including frequency- 
dependent, lossy electric permittivity and magnetic permeability tensors; 

• objects separated by vacuum or by a medium with uniform, frequency-dependent 
isotropic permittivity and permeability; 

• zero or nonzero temperature; 

• and objects outside of one another or enclosed in each other. 

These ideas build on a range of previous related work, an inevitably incomplete 
subset of which is briefly reviewed here: Scattering theory methods were first ap- 
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plied to the parallel plate geometry, when Kats reformulated Lifshitz theory in terms 
of reflection coefficients [50|. Jaekel and Reynaud derived the Lifshitz formula us- 
ing reflection coefficients for lossless infinite plates [48 1 and Genet, Lambrecht, and 
Reynaud extended this analysis to the lossy case [ 38 1 . Lambrecht, Maia Neto, and 
Reynaud generalized these results to include non-specular reflection |58 1. 

Around the same time as Kats's work, Balian and Duplantier developed a mul- 
tiple scattering approach to the Casimir energy for perfect metal objects and used 
it to compute the Casimir energy at asymptotically large separations (4J [3) at both 
zero and nonzero temperature. In their approach, information about the conductors 
is encoded in a local surface scattering kernel, whose relation to more conventional 
scattering formalisms is not transparent, and their approach was not pursued further 
at the time. One can find multiple scattering formulas in an even earlier article by 
Renne [81 1, but scattering is not explicitly mentioned, and the technique is only used 
to rederive older results. 

Another scattering-based approach has been to express the Casimir energy as an 
integral over the density of states of the fluctuating field, using the Krein formula 
ll56l l57l |6) to relate the density of states to the S-matrix for scattering from the 
ensemble of objects. This S-matrix is difficult to compute in general. In studying 
many-body scattering, Henseler and Wirzba connected the 5-matrix of a collection 
of spheres [47 1 or disks |93 1 to the objects' individual 5-matrices, which are easy to 
find. Bulgac, Magierski, and Wirzba combined this result with the Krein formula to 
investigate the scalar and fermionic Casimir effect for disks and spheres lfTTl[T0ll94ll . 
Casimir energies of solitons in renormalizable quantum field theories have been 
computed using scattering theory techniques that combine analytic and numerical 
methods ll44l . 

Bordag, Robaschik, Scharnhorst, and Wieczorek [7, 82 1 introduced path integral 
methods to the study of Casimir effects and used them to investigate the electro- 
magnetic Casimir effect for two parallel perfect metal plates. Li and Kardar used 
similar methods to study the scalar thermal Casimir effect for Dirichlet, Neumann, 
and mixed boundary conditions |62, 63 1. The quantum extension was developed fur- 
ther by Golestanian and Kardar l4T1l42l and was subsequently applied to the quan- 
tum electromagnetic Casimir effect by Emig, Hanke, Golestanian, and Kardar, who 
studied the Casimir interaction between plates with roughness [ 3 1 1 and between de- 
formed plates [32]. (Techniques developed to study the scalar Casimir effect can be 
applied to the electromagnetic case for perfect metals with translation symmetry in 
one spatial direction, since then the electromagnetic problem decomposes into two 
scalar ones.) Finally, the path integral approach was connected to scattering theory 
by Emig and Buescher lfL2ll . 

Closely related to the work we present here is that of Kenneth and Klich, who 
expressed the data required to characterize Casimir fluctuations in terms of the tran- 
sition T-operator for scattering of the fluctuating field from the objects BTI . Their 
abstract representation made it possible to prove general properties of the sign of 
the Casimir force. In Refs. [28, 29|, we developed a framework in which this ab- 
stract result can be applied to concrete calculations. In this approach, the T-operator 
is related to the scattering amplitude for each object individually, which in turn is 
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expressed in an appropriate basis of multipoles. While the T-operator is in gen- 
eral "off-shell," meaning it has matrix elements between states with different spatial 
frequencies, the scattering amplitudes are the "on-shell" matrix elements of this op- 
erator between states of equal spatial frequency!*] So, it is not the T-operator itself 
that connects, say, outgoing and standing waves in the case of outside scattering but 
its on-shell matrix elements, the scattering amplitudes. In this approach, the objects 
can have any shape or material properties, as long as the scattering amplitude can 
be computed in a multipole expansion (or measured). The approach can be regarded 
as a concrete implementation of the proposal emphasized by Schwinger |90| that 
the fluctuations of the electromagnetic field can be traced back to charge and cur- 
rent fluctuations on the objects. This formalism has been applied and extended in a 
number of Casimir calculations l52l l67l l68l l80l l40l l9ll . 

The basis in which the scattering amplitude for each object is supplied is typ- 
ically associated with a coordinate system appropriate to the object. Of course a 
plane, a cylinder, or a sphere would be described in Cartesian, cylindrical, or spher- 
ical coordinates, respectively. However, any compact object can be described, for 
example, in spherical coordinates, provided that the matrix of scattering amplitudes 
can be either calculated or measured in that coordinate system. There are a limited 
number of coordinate systems in which such a partial wave expansion is possible, 
namely those for which the vector Helmholtz equation is separable. The translation 
matrices for common separable coordinate systems, obtained from the free Green's 
function, are supplied in Appendix C of reference [75 1. For typical cases, the final 
computation of the Casimir energy can be performed on a desktop computer for a 
wide range of separations. Asymptotic results at large separation can be obtained 
analytically. 

The primary limitation of the method is on the distance between objects, since the 
basis appropriate to a given object may become impractical as two objects approach. 
For small separations, sufficient accuracy can only be obtained if the calculation is 
taken to very high partial wave order. (Vastly different scales are problematic for nu- 
merical evaluations in general.) In the case of two spheres, the scattering amplitude 
is available in a spherical basis, but as the two spheres approach, the Casimir en- 
ergy is dominated by waves near the point of closest approach |89|. As the spheres 
come into contact an infinite number of spherical waves are needed to capture the 
dominant contribution. A particular basis may also be fundamentally inappropriate 
at small separations. For instance, if the interaction of two elliptic cylinders is ex- 
pressed in an ordinary cylindrical basis, when the elliptic cylinders are close enough, 
the smallest circular cylinder enclosing one may not lie outside the smallest circular 
cylinder enclosing the other. In that case the cylindrical basis would not "resolve" 
the two objects (although an elliptic cylindrical basis would). Finally, for a variety 
of conceptual and computational reasons, we are limited to linear electromagnetic 
response. 



1 Because of this relationship, these scattering amplitudes are also referred to as elements of the 
T -matrix. In standard conventions, however, the T-matrix differs from the matrix elements of the 
T-operator by a basis-dependent constant, so we will use the term "scattering amplitude" to avoid 
confusion. 
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In spirit and in mathematical form our final result resembles similar expressions 
obtained in surface integral equation methods used in computational electrodynam- 
ics l2ll . Using such a formulation, in which the unknowns are currents and fields 
on the objects, one can compute the Casimir energy using more general basis func- 
tions, e.g., localized basis functions associated with a grid or mesh, giving rise to 
finite element and boundary elements methods [80|. 

In addition to an efficient computational approach, the scattering formalism has 
provided the basis for proving general theorems regarding Casimir forces. The 
seemingly natural question whether the force is attractive or repulsive turns out to 
be an ill-defined or, at least, a tricky one on closer inspection. When, for example, 
many bodies are considered, the direction of the force on any one object depends, of 
course, on which other object's perspective one takes. Even for two objects, "attrac- 
tive" forces can be arranged to appear as a "repulsive" force, as in the case of two 
interlocking combs [84]. To avoid such ambiguous situations one can restrict oneself 
to analyzing two objects that are separable by a plane. Even here, it has turned out 
that a simple criterion for the direction of the force could not be found. Based on var- 
ious calculations for simple geometries it was thought that the direction of the force 
can be predicted based on the relative permittivities and permeabilities of the objects 
and the medium. Separating materials into two groups, with (i) permittivity higher 
than the medium or permeability lower than the medium (e > Em and ji < jXm), or 
(ii) the other way around (e < Em and ji > jj,m), Casimir forces had been found to be 
attractive between members of the same group and repulsive for different types in 
the geometries considered. However, a recent counterexample [61 1 shows that this 
is not always true. A rigorous theorem, which states that Casimir forces are always 
attractive, exists only for the special case of mirror symmetric arrangements of ob- 
jects. It was proven first with a T-operator formalism [51], similar to our approach 
used here, and later using reflection positivity (3). We have taken an alternative 
characterization of the force to be fundamental, namely, whether it can produce a 
stable equilibrium [77]. Here, the categorization of materials into the two groups is 
meaningful since objects made of materials of the same type cannot produce stable 
levitation. One practical consequence of this theorem is that it reveals that many 
current proposals for producing levitation using metamaterials cannot succeed. 

To illustrate the general formulation, we provide some sample applications. We 
include an analysis of the forces between two cylinders or wires [76 1 and a cylinder 
and a plate lf33l l76l 1751 . The Casimir interaction of three bodies is presented sub- 
sequently; it reveals interesting multibody effects ll78l l76l l85l . The Casimir torque 
of two spheroids is discussed as well [30]. Furthermore, we analyze the Casimir 
effect for a parabolic cylinder opposite a plate when both represent perfect metal 
material boundary conditions [45 1. We find that the Casimir force does not vanish 
in the limit of an infinitesimally thin parabola, where a half plate is arranged above 
an infinite plate, and we compute the edge effect. Another type of geometries that is 
treated here consists of a finite sphere or a small spheroid inside a spherical metallic 
cavity l95l . 

This chapter is organized as follows: First, we sketch the derivation of the 
Casimir interaction energy formula Eq. ( 39 1 in Section [2] Next, the theorem re- 
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garding stability is derived in Section [3] Finally, in Section [4] sample applications 
are presented. 



2 General theory for Casimir interactions 

This section has been adapted from a longer article, reference ll75l . Many tech- 
nical details and extensive appendices have been omitted to fit the format of this 
book. 



2.1 Path integral quantization 

2.1.1 Electromagnetic Lagrangian 

We consider the Casimir effect for objects without free charges and currents but 
with nonzero electric and magnetic susceptibilities. The macroscopic electromag- 
netic Lagrangian density is 

J*f=i(E D-B H). (1) 

The electric field E(f,x) and the magnetic field B(?,x) are related to the funda- 
mental four- vector potential A^ by E = — c~ l d t A — VA° and B = V x A. We treat 
stationary objects whose responses to the electric and magnetic fields are linear. For 
such materials, the D and B fields are related to the E and H fields by the convo- 
lutions D(f,x) = fZ K dt f e(t f ,x)E(t — t f ,x) and B(f,x) = fZ K dt'n(t f ,x)TL(t — t f ,x) 
in time, where e(t',x) and ju(f',x) vanish for t' < 0. We consider local, isotropic 
permittivity and permeability, although our derivation can be adapted to apply to 
non-local and non-isotropic media simply by substituting the appropriate non-local 
and tensor permittivity and permeability functions. A more formal derivation of our 
starting point Eq. (fill, which elucidates the causality properties of the permeability 
and permittivity response functions, is given in Appendix A of reference [75 1. 

We define the quantum-mechanical energy through the path integral, which sums 
all configurations of the electromagnetic fields constrained by periodic boundary 
conditions in time between and T. Outside of this time interval the fields are 
periodically continued. Substituting the Fourier expansions of the form E(f,x) = 
^™ = _ 00 E(o)„,x)e _ " B "' with (O n = 2nn/T, we obtain the action 

S(T) = - J dt fdx(E-D-B-U) = -T £ /*dx(E*-£E-B*-jU _1 B), (2) 



where e, E, /i, and B on the right-hand side are functions of position x and frequency 

!_ 



(O n , and we have used D(a),x) = e(a),x)E(a),x) and H((»,x) = ,^, 8(0), x) 
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From the definition of the fields E and B in terms of the vector potential A^, we 
have VxE = /— B, which enables us to eliminate B in the action, 



1 



S(T) = -T £ hlx 



E* 



to, 



Vx V; 



E 



r E*-VE 



(o- 



(3) 



where 



a, 



I^-(l-e(»„,x)) + Vx 



1 



1 Vx 



(4) 



k jU(C0„,X 

is the potential operator and we have restored the explicit frequency dependence of 
£ and fi. The potential operator is nonzero only at those points in space where the 
objects are located (e ^ 1 or jj. ^ 1). 

In the functional integral we will sum over configurations of the field A^ . This 
sum must be restricted by a choice of gauge, so that it does not include the infinitely 
redundant gauge orbits. We choose to work in the gauge A = 0, although of course 
no physical results depend on this choice. 



2.1.2 Casimir energy from Euclidean action 



We use standard tools to obtain a functional integral expression for the ground state 
energy of a quantum field in a fixed background described by V(Ct),x). The overlap 
between the initial state |E„) of a system with the state |E&) after time T can be ex- 
pressed as a functional integral with the fields fixed at the temporal boundaries [36], 



(E 6 | e -' H '*|E a )= / ^A| E( , =0)=Efl ^ 

E(,=T)=E b 



ism 



(5) 



where S(T) is the action of Eq. (pi with the time integrals taken between zero and 
T, and H is the corresponding Hamiltonian. 

If the initial and final states are set equal and summed over, the resulting func- 
tional integration defines the Minkowski space functional integral 



2?{T) =Y,(E a \e- iHT/n K) =Xxe- im l % = [ @Ae 



AS(T) 



(6) 



which depends on the time T and the background potential V(ffl,x). The partition 
function that describes this system at temperature 1 /)3 is defined by 



Z(j8) = ^(-mj8)=tre-^ H , 
and the free energy F of the field is 

F(j8) = -ilogZ(j8). 



(7) 



(8) 
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The limit j3 — > °° projects the ground state energy out of the trace, 

^ =n/3=°°) = -jim,!rlogZ(j3). (9) 

The unrenormalized energy £q generally depends on an ultraviolet cutoff, but cutoff- 
dependent contributions arise from the objects individually l43ll44ll and do not de- 
pend on their separations or orientations. Such terms can remain after ordinary QED 
renormalization if objects are assumed to constrain electromagnetic waves with ar- 
bitrarily high frequencies (for example, if the fields are forced to vanish on a sur- 
face). Such boundary conditions should be regarded as artificial idealizations; in 
reality, when the wavelengths of the electromagnetic waves become shorter than the 
length scales that characterize the interactions of the material, the influence of the 
material on the waves vanishes |43l . Accordingly, the potential V should vanish for 
real materials in the high-frequency limit. In any event these cutoff dependences are 
independent of the separation and orientation of the objects, and since we are only 
interested in energy differences, we can remove them by subtracting the ground 
state energy of the system when the objects are in some reference configuration. 
In most cases we take this configuration to be when the objects are infinitely far 
apart, but when calculating Casimir energies for one object inside another, some 
other configuration must be used. We denote the partition function for this reference 
configuration by Z. In this way we obtain the Casimir energy, 

<f = -lim-LogZ(j3)/Z(j3). (10) 

Throughout our calculation of S, we will thus be able to neglect any overall factors 
that are independent of the relative positions and orientations of the objects. 

By replacing the time T by —infi, we transform the Minkowski space functional 
integral 3f(T) into the partition function Z(j3). In A = gauge, the result is simply 
to replace the frequencies co„ = ^ in Eq. (|4k by OSr = icK n , where K n is the n th 

Matsubara frequency divided by c. (In other gauges the temporal component A of 
the vector field must be rotated too.) 

The Lagrangian is quadratic, so the modes with different K„ decouple and the 
partition function decomposes into a product of partition functions for each mode. 
Since the electromagnetic field is real, we have E*(o) = E(— (0) on the real axis. 
We can thus further simplify this decomposition on the imaginary axis by consid- 
ering K > only, but allowing E and E* to vary independently in the path integral. 
Restricting to positive K is possible because the response functions e(icK,x) and 
jj,(icK,x) are invariant under a change of sign in icK, as shown in Appendix A of 
Ref. [75]. In the limit j8 —¥ °°, the sum £n>0 turns into an integral -^ JJ° dK, and 

we have 

he C°° 

4 = -—/ dK\ogZ(K), (11) 

where 
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Z(k)= J 0A0A*exp 



-J3 / dxE* ■ ( I+— Vx Vx ]E + — E*-V(zcJC,x)E 



(12) 



V(/cK-,x)=lK- z (e(/cK-,x)-l) + Vx — --1 Vx. (13) 

\H(lCK,X) ) 

The potential Y(icK,x) is real for real K, even though e and jj. can have imaginary 



parts for real frequencies CO. Our goal is now to manipulate Z(k) in Eq. ( 12 1 so that 
it is computable from the scattering properties of the objects. 



2.2 Green's function expansions and translation formulas 

The free Green's function and its representations in various coordinate systems are 
crucial to our formalism. The free electromagnetic field (V = 0) obeys equations of 
motion obtained by extremizing the corresponding action, Eq. p}, 



CO 2 



c 



2 



VxVx E(o),x)=0. (14) 



We will employ the electromagnetic dyadic Green's function Go, defined by 

, rVxVx]Gn(a),x,x') = I<5 (3) (x-x'), (15) 



CO 2 



written here in the position space representation. The Green's function has to be 
the retarded one, not only on physical grounds, but also as a consequence of the 
imaginary-frequency formalism, just as is the case for the response functions e and 
jj,. It is the retarded response functions that are analytically continued in the fre- 
quency domain to positive imaginary frequency, as shown in Appendix A of refer- 
ence IT751 . 

The representation of the free Green's function, which we need, employs the 



"regular" and "outgoing" solutions to the differential equation, Eq. ( 14 1, 

E^(co,x) = (x\E^(co)), ET(co,x) = (x\ET(co)), (16) 

represented formally by the eigenstate kets |E„ g (a))) and |Ea S (a))), where the gen- 
eralized index a labels the scattering channel, including the polarization. For ex- 
ample, for spherical wave functions it represents the angular momentum quantum 
numbers (l,m) and the polarization E or M. There are six coordinate systems in 



which the vector wave equation ( 14 1 can be solved by separation of variables and 
vector wave functions appropriate to that coordinate system can be constructed [70]. 
The labels "regular" and "outgoing" denote, respectively, the wave functions' non- 
singular behavior at the origin or 'outward' direction of energy transport along one 
of the coordinate system's axes. Let us call the coordinate, along which the latter 
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wave functions are outgoing %\ and the other coordinates ^2 and §3. We will usually 
work on the imaginary co-axis, in which case we will encounter the corresponding 
modified special functions. 

The free Green's function can be expanded in tensor products of these wave 
functions, 



)(ffl,x,x')=£C a i 



to 



E^(©,Si,&,&)®ES g >,&&3) if ^i(x) > Si?) 

K % {<0^2&)®K*M V ^) if^(x) <^'(X') ; 

(17) 
E™ is the same as E° ut except the functional dependence on %\ is complex conju- 
gated, making the wave function 'incoming'. A list of Green's function expansions 
in various common bases, including the normalization constant, C a (co), is given in 
Appendix B of Ref. [75 J. The wave functions that appear in the series expansion 
of the free Green's functions in Eq. ( fl7| satisfy wave equations with frequency CO. 
As we will see in Sec. |2.3| the ability to express the Casimir energy entirely in 
terms of an "on-shell" partial wave expansion with fixed CO will greatly simplify our 
calculations. 

We will also use the free Green's function in another representation to combine 
the scattering amplitudes for two different objects. In this calculation the one argu- 
ment of the Green's function will be located on each object. As long as the pair of 
objects can be separated in one of the separable coordinate systems by the surface 
4i = S =const., we can distinguish an inside object which lies entirely inside the 
surface (<jji < S) and an outside object (<jji > 3), see Fig. [2] Then, we can expand 





Fig. 2 Geometry of the outside (left) and inside (right) configurations. The dotted lines show 
surfaces separating the objects on which the radial variable is constant. The translation vector 
X, ; = x, — Xj = — Xy; describes the relative positions of the two origins. 
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the free Green's function, when one argument, say x, lies on object i and the other 
argument, say x', lies on object j, we expand Go(/cK",x,x') in terms of coordinates 
x, and x' that describe each point relative to the origin of the body on which it lies. 
Which of the following expansions is appropriate for a particular pair of objects de- 
pends whether objects i and j are outside of one another, or one object is inside the 
other, 

icK,x,x') — 

' E^ s (k-,x,) <g> <%K (fc)E 1 „ eg *(K-,x^) if i and j are outside each other 



EW 

a.p 



-^vee, n jAi , lr ,i„/ , x [if / is inside /', or 

E r „ cg (K,Xi )®VA {k)WS* (k,x'A I ...... '.,. . ■ x 

v ' «P V ' P v J J I if i is below j (plane wave basis) 

ET(K^i) ® «(«>,£) ( ^ iS |" S . ide (' * , . . 

" v ' a P v ' P v •'•' ^ if j is below i (plane wave basis) 

(18) 



where Wi' a (K) = Y^' R (K) ^ a )1[ and C n is the normalization constant defined in 



»*«S(*) = ^5V)j$g 



Eq. ( 17 I, The expansion can be written more compactly as 



/ (E reg (K-)l\ 
Go(icv) = E(-W) (|E« g «) |Er(v)))X^(!c) f y n( ^' j , (19) 

where the X matrix is defined, for convenience, as the negative of the matrix con- 
taining the translation matrices, 



In Eq. ( 19 1 the bras and kets are to be evaluated in position space in the appropriately 
restricted domains and only one of the three submatrices is nonzero for any pair of 
objects i and j as given in Eq. (|18), The translation matrices for various geometries 



are provided in Appendix C of reference 11751 . 



2.3 Classical scattering of electromagnetic fields 

In this section, we summarize the key results from scattering theory needed to com- 
pute the scattering amplitude of each body individually. In the subsequent section 
we will then combine these results with the translation matrices of the previous 
section to compute Z(k). 

By combining the frequency-dependent Maxwell equations, one obtains the vec- 
tor wave equation 

0) 2 

(Ho + V(©,x))E(a>,x) = -^-E(o),x), (21) 
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where 

H = Vx Vx, 

co 2 ( 1 \ (22) 

which is the same potential operator as the one obtained by rearranging the La- 
grangian (see Eq. Q). 

The Lippmann-Schwinger equation 

|E) = |Eo)-G V|E) (23) 



sor Green's function discussed in Sec 



expresses the general solution to Eq. (|2T|). Here, Go is the free electromagnetic ten- 

Eo) obeys 



2.2 and the homogeneous solution 



f-^I+Ho) |Eo) = 0. We can iteratively substitute for |E) in Eq. (23 1 to obtain 
the formal expansion 



|E) = |E ) - GoV|Eo) + GoVG V|E) - . . . 
= |E )-GoT|Eo), 

where the electromagnetic T-operator is defined as 

T= VfT I^=VGG i, (25) 

and G is the Green's function of the full Hamiltonian, (- ^I + H + v] G = I. We 
note that T, Go, and G are all functions of frequency CO and non-local in space. 



As can be seen from expanding T in Eq. (25 i in a power series, T(e),x,x') = 
(x|T(ffl)|x / ) is zero whenever x or x' are not located on an object, i.e., where V(o),x) 
is zero. This result does not, however, apply to 

T _1 =G +V _1 , (26) 

because the free Green's function is nonlocal. The potential V(a),x) which appears 



inEq. ( 22 lis the coordinate space matrix element of V, (x|V|x') = V(fi),x)5(x— x'), 
which can be generalized to the case where V is non-local, (x|V|x') = V(o,x,x'). 
Note that whether V is local or non-local, its matrix elements vanish if x and x' are 
on different objects or if either x or x' is outside of the objects. The definition of V~ 
is natural, (x|V _1 |x') = V _1 (a>,x)5(x — x') (and similarly for the non-local case) 
when x and x' are on a single object, which is the only case that enters our analysis. 

Next we connect the matrix elements of the T-operator between states with equal 
CO to the scattering amplitude J?. In our formalism, only this restricted subset of T- 
operator matrix elements is needed in the computation of the Casimir energy. 

By the choice of the homogeneous solution, |Eo), is regular or outgoing, we can 
distinguish two physically different processes. In the former case, the object scatters 
the regular wave outward and modifies the amplitude of the imposed regular wave 
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Fig. 3 The scattering waves for outside scattering (left panel) and inside scattering (right panel). 
In both cases the homogeneous solution Eo(co) is shown in bold. For outside scattering, the ho- 
mogeneous solution is a regular wave, which produces a regular wave inside the object and an 
outgoing wave outside the object. For inside scattering, the homogeneous solution is an outgoing 
wave, which produces a regular wave inside the object and an outgoing wave outside the object. 



functions inside, a situation we refer to as outside scattering (left panel of Fig.[3||^Jln 
the latter case, the object modifies the amplitude of the transmitted wave and partly 
reflects it as a regular wave inside (inside scattering, right panel of Fig. [3]). 'Outside' 
and 'inside' are distinguished by surfaces ^i=constant, as before. Here we treat the 
outside scattering case, and refer the reader to Ref. [75] the inside case and further 
details. The expansion in Eq. ( fTT) allows us to express Eq. ( [24| as 



E(a),x)=E r a eg (a),x)-£E£ ut (c0,x) 



(27) 



x / C J3 (fi))E™ g *(©,x') ■T(m,x',x")E% s {a),x")dx'dx". 

at points x outside a surface £,i =constant enclosing the object. The equation can be 
written in Dirac notation, again with the condition that the domain of the functional 
Hilbert space is chosen appropriately to the type of solution, 

|E(o))) = |E r a eg (a)))+£|E£ ut (a))) x (-l)C p (co)(E^(co)\nco)\K"(co)), (28) 
J3 " v ' 

which defines <^a e a as the exterior/exterior scattering amplitude (the one evaluated 
between two regular solutions). We will use analogous notation in the other cases 
below. 

At coordinates x "far enough inside" a cavity of the object, meaning that x has 
smaller £,\ than any point on the object, the field E is given by 



Alternatively, we can set up asymptotically incoming and outgoing waves on the outside and 
regular waves inside. The amplitudes of the outgoing waves are then given by the 5-matrix, which 
is related to the scattering amplitude .$- by & = (S — I)/2. Although these two matrices carry 
equivalent information, the scattering amplitude will be more convenient for our calculation. 
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|E(cd)) = \E^(co)) + £|E re ») x (-l)Cp(co)(E^(co)\T(co)\E^(co)), (29) 
p s * ' 

where again the free states are only defined over the appropriate domain in position 
space, and jF' e indicates the interior/exterior scattering amplitude. 

We have obtained the scattering amplitude in the basis of free solutions with fixed 
CO. Since one is normally interested in the scattering of waves outside the object, the 
scattering amplitude usually refers to JF ee . We will use a more general definition, 
which encompasses all possible combinations of inside and outside. The scattering 
amplitude is always "on-shell," because the frequencies of both the operator and the 
states is CO. As a result, it is a special case of the T-operator, which can connect wave 
functions with different CO. 

We find it convenient to assemble the scattering amplitudes for inside and outside 
into a single matrix, 

' E^(K)\T(icK)\E^(K)) (K*(K)\T(icK)\E°«(K))\ 

[E%(K)\T(icK)\Ef(K)) (E%(K)\T(icK)\E°«(K)) y 

(30) 



¥(K) = 

= (-l)C a (K) 



where we have set CO — icK, since this is the case we use. For simplificity we define 



2.4 Casimir free energy in terms of the scattering amplitudes 

With the tools of the previous two sections, we are now able to re-express the 



Euclidean electromagnetic partition function of Eq. ( 12 1 in terms of the scattering 



theory results derived in Section 2.3 for imaginary frequency. We exchange the fluc- 
tuating field A, which is subject to the potential V(icK, x), for a free field A', together 
with fluctuating currents J and charges — ~ V • J that are confined to the objects 3 



We multiply and divide the partition function Eq. ( fT2| by 



W= /^J^r| obj exp 



-J3 /dxJ*(x)-V _1 (jcK-,x)J(x 



= detV(icK,x,x'), 

(31) 



where | ob j indicates that the currents are defined only over the objects, i.e. the do- 
main where V is nonzero and therefore V~' exists. 



3 The sequence of two changes of variables is known as Hubbard-Stratonovich transformation in 
condensed matter physics. 
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We then change variables in the integration, J(x) = J'(x) + ^V(/cfC,x)E(x), and 
a second time, E(jcJC,x) = E'(/cfC,x) — iK J dx' Go(icK,x,x')J (x 1 ) and analogously 
for J* and E*, to obtain 



zM-l/sWU, 



exp 



-P / dxdx'j'*(x)-(G (icK,x,x')+Y- l {icK,x,x'))j'(x') 



(32) 



where 



Z = /0A'0A'*exp 



-P / rfxE'*(x) 



r Vx Vx jE'(x) 



(33) 



is the partition function of the free field, which is independent of the objects. In 



Z(k), current fluctuations replace the field fluctuations of Eq. ( 12 1. The interaction 



of current fluctuations on different objects is described by the free Green's func- 
tion Go(icK,x,x') alone. The inverse potential penalizes current fluctuations if the 
potential is small. 

To put the partition function into a suitable form for practical computations, we 
use the results of the previous sections to re-express the microscopic current fluc- 
tuations as macroscopic multipole fluctuations, which then can be connected to the 
individual objects' scattering amplitudes. This transformation comes about naturally 
once the current fluctuations are decomposed according to the objects on which they 
occur and the appropriate expansions of the Green's function are introduced. We 
begin this process by noticing that the operator in the exponent of the integrand in 



Eq. ( 32 1 is the negative of the inverse of the T-operator (see Eq. (26 1), and hence 



Z(K-)=Z detV _1 (/cK-,x,x')detT(/cK-,x,x') 



(34) 



which is in agreement with a more formal calculation: Since Zq = detGo(/cfC,x,x') 
and Z(k) = detG(/cfC,x,x'), we only need to take the determinant of Eq. (25 1 to 
arrive at the result of Eq. (J34]i. 

Both Zo and detV _1 (?cK",x) are independent of the separation of the objects, 
since the former is simply the free Green's function, while the latter is diagonal 
in x. Even a nonlocal potential V(icJC,x,x') only connects points within the same 
object, so its determinant is also independent of the objects' separation. Because 
these determinants do not depend on separation, they are canceled by a reference 
partition function in the final result. We are thus left with the task of computing the 
determinant of the T-operator. 



As has been discussed in Sec. 2.3 the T-operator T(z'cfC,x,x') is not diagonal in 



the spatial coordinates. Its determinant needs to be taken over the spatial indices x 
and x', which are restricted to the objects because the fluctuating currents J(x) in 
the functional integrals are zero away from the objects. This determinant also runs 
over the ordinary vector components of the electromagnetic T operator. 
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A change of basis to momentum space does not help in computing the determi- 
nant of the T-operator, even though it does help in finding the determinant of the 
free Green's function, for example. One reason is that the momentum basis is not 
orthogonal over the domain of the indices x and x', which is restricted to the objects. 
In addition, a complete momentum basis includes not only all directions of the mo- 
mentum vector, but also all magnitudes of the momenta. So, in the matrix element 
(Ek|T(o)|E[;/) the wave numbers k and k' would not have to match, and could also 
differ from co/c. That is, the matrix elements could be "off-shell." Therefore, the 
T-operator could not simply be treated as if it was the scattering amplitude, which 
is the on-shell representation of the operator in the subbasis of frequency CO (see 
Sec. |2.3| >, and is significantly easier to calculate. Nonetheless, we will see that it is 
possible to express the Casimir energy in terms of the on-shell operator only, by 
remaining in the position basis. 



From Eq. (25 I, we know that the inverse of the T-operator equals the sum of the 
free Green's function and the inverse of the potential. Since the determinant of the 
inverse operator is the reciprocal of the determinant, it is expedient to start with the 
inverse T-operator. We then separate the basis involving all the objects into blocks 
for the n objects. In a schematic notation, we have 



/[(xilTrVj)] 



[(x|T-V>]= [telGoft) 



V" 



[<Xi|Go|*£ 



[(x 2 |T 2 -^)] 



7 



(35) 



where the // h submatrix refers to x G object i and x' £ object j and x,- represents a 
point in object i measured with respect to some fixed coordinate system. Unlike the 



position vectors in Sec. 2.2 at this point the subscript of x, does not indicate the ori- 
gin with respect to which the vector is measured, but rather the object on which the 
point lies. Square brackets are used to remind us that we are considering the entire 
matrix or submatrix and not a single matrix element. We note that the operators T 
and Go are functions of icK, but for simplicity we suppress this argument throughout 
this derivation. When the two spatial indices lie on different objects, only the free 
Green's function remains in the off-diagonal submatrices, because (x;|V _1 |x'-) = 
for i^j. 

Next, we multiply T _1 by a reference T-operator Too without off-diagonal sub- 
matrices, which can be interpreted as the T-operator at infinite separation, 



(x|TooT 



-l 



x")l 



[<xi| 



[;</x 2 (x 2 |T 2 |x 2 )(x 2 |G |x'/)] 



[J^ 1 (x 1 |T 1 |x;)(x;|Go|x 2 Q 



(x 2 |x' 2 ' 



7 



(36) 



Each off-diagonal submatrix [/rfx|-(x,-|T,-|x|-)(x;|Go|x'')] is the product of the T- 
operator of object i, evaluated at two points x* and x'- on that object, multiplied 
by the free Green's function, which connects xj to some point x" on object /'. 
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Now we shift all variables to the coordinate systems of the objects on which 
they lie. As a result, the index on a position vector x,- now refers to the object i 
on which the point lies and to the coordinate system with origin 0\ in which the 



vector is represented, in agreement with the notation of Sec. 2.2 The off-diagonal 
submatrices in Eq. ( 36 1 can then be rewritten using Eq. (fT9]i as, 



((x,|T,K eg M> (xi|T,-|Er( J c)))X^ 




(-q,(K)) 



(37) 



The matrix [(x|TooT ^x")] has the structure I + AB. Using Sylvester's determi- 
nant formula det(I + AB) = det(I + BA), we see that the determinant is unchanged 



if we replace the off-diagonal submatrices in Eq. ( 36 1 by 



K-i)Q(k) 



'^{k)\T^J(k)) (E™ 8 ( K )|T,-|E^(«:)> \ y 
(E£(K)|T,|E r ; s (K)) (ES(K)|T,-|E^(f)> I ?'■ 



(38) 



With this change, the diagonal submatrices in Eq. ( |36| ) become diagonal in the par- 
tial wave indices rather than in position space. The matrix elements of the T-operator 
are the scattering amplitudes, which can be obtained from ordinary scattering cal- 



culations, as demonstrated in Sec. 2.3 The first matrix in Eq. (38 1, including the 
prefactor (— l)C a (fc), is F,-(jc), the modified scattering amplitude of object i, de- 
fined in Eq. ([30}. 



(39) 



Putting together Eqs. ( fTT) , (12) , ffify , and ( |36| , we obtain 



£■■ 



he 
2nJo 



djdogdet( 



')■ 



where 



?-i 



-12 



X 21 F^ 1 



[ 13 

-23 



(40) 



and M«, is a block diagonal matrix diag(Fi 
Using the block determinant identity 



det 



det(A)det( 



.-It 



:det(D)det( 



1,-1 



C) , (41) 



we can simplify this expression for the case of the interaction between two objects, 



§■■ 



he 
2nJo 



dKlogdet(l-F a X ab F h l 



•ho. 



(42) 



Usually, not all of the submatrices of F and X are actually needed for a computa- 
tion. For example, if all objects are outside of one another, only the submatrices ^ ee 
of the scattering amplitude that describe outside reflection are needed. If there are 
only two objects, one inside another, then only the inside reflection submatrix jF" 
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of the outside object and the outside reflection submatrix & ee of the inside object 
are needed. 

In order to obtain the free energy at nonzero temperature instead of the ground 
state energy, we do not take the limit /3 — >• °° in Eq. d9b. Instead, the integral || /J° dK 
is replaced everywhere by A £j,, where ck„ = j£- with n = 0, 1,2, 3 . . . is the nth 
Matsubara frequency. A careful analysis of the derivation shows that the zero fre- 
quency mode is weighted by 1/2 compared to the rest of the terms in the sum; this 
modification of the sum is denoted by a prime on the summation symbol. The fac- 
tor of 1/2 comes about because the fluctuating charges or currents have to be real 
for zero frequency. Thus, for Kq, the expressions on the right hand side of Eq. ( [34} 
should be placed under a square root. (For a complex field, both signs of the integer 
n would be included separately, and n = would be included once, with the normal 
weight.) 

If the medium between the objects is not vacuum but instead has permittivity 
EatiicK) and magnetic permeability plm{wk) different from unity, then the free 
Green's function is multiplied by pLm{ick), and its argument K is replaced by 
nM{icK)K, where tiM{icK) = WEm (ick) \Im (ick) is the medium's index of refrac- 
tion. Effectively, this change just scales all frequency dependencies in the transla- 
tion matrices X(jc), which become X(riM{icK)K). Furthermore, the scattering am- 
plitudes absorb the factor jj.m(ick) from the free Green's function and change non- 
trivially, i.e. not just by some overall factor or a scaling of the frequency. They 
have to be computed with the nonzero electric and magnetic susceptibilities of the 
medium. 



3 Constraints on stable equilibria 



Before presenting particular applications of the Casimir energy expression in Eq. ( 39 1 
we consider some general properties of electrodynamic Casimir interactions here. 
This section has been adapted from a letter, which is co-authored by two of us f77l . 

As described in the Introduction, some general statements about the attractive or 
repulsive nature of Casimir forces can be made on the basis of the relative permittiv- 
ity and permeability of objects and the medium they are immersed in. But the sign 
of the force is largely a matter of perspective, since attractive forces can be easily 
arranged to produce repulsion along a specific direction, e.g., as in Ref. |84|. In- 
stead, we focus on the question of stability, see Fig. [4] which is more relevant to the 
design and development of MEMs and levitating devices. We find that interactions 
between objects within the same class of material (as defined in the Introduction) 
cannot produce stable configurations. 

Let us take a step back and consider the question of stability of mechanical equi- 
libria in the realm of electromagnetism. Earnshaw's theorem [26 1 states that a col- 
lection of charges cannot be held in stable equilibrium solely by electrostatic forces. 
The charges can attract or repel, but cannot be stably levitated. While the stability 
of matter (due to quantum phenomena) is a vivid reminder of the caveats to this 
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R 




Fig. 4 The Casimir energy is considered for objects with electric permittivity e,(tt),x) and mag- 
netic permeability /i,(fi),x), embedded in a medium with uniform, isotropic, £m((q) and jXm{w). 
To study the stability of object A, the rest of the objects are grouped in the combined entity R. The 
stability of the position of object A is probed by displacing it infinitesimally by vector d. 



theorem, it remains a powerful indicator of the constraints to stability in electrostat- 
ics. An extension of Earnshaw's theorem to polarizable objects by Braunbek [9, 8| 
establishes that dielectric and paramagnetic (e > 1 and jj. > 1) matter cannot be sta- 
bly levitated by electrostatic forces, while diamagnetic (fi < 1) matter can. This is 
impressively demonstrated by superconductors and frogs that fly freely above mag- 
nets [37]. If the enveloping medium is not vacuum, the criteria for stability are mod- 
ified by substituting the static electric permittivity Em and magnetic permeability jj,m 
of the medium in place of the vacuum value of 1 in the respective inequalities. In 
fact, if the medium itself has a dielectric constant higher than the objects (e < Em), 
stable levitation is possible, as demonstrated for bubbles in liquids (see Ref. ||49l , 
and references therein). For dynamic fields the restrictions of electrostatics do not 
apply; for example, lasers can lift and hold dielectric beads with index of refraction 
n = y/Ejl > 1 |0Q] . In addition to the force which keeps the bead in the center of the 
laser beam there is radiation pressure which pushes the bead along the direction of 
the Poynting vector. Ashkin and Gordon have proved that no arrangement of lasers 
can stably levitate an object just based on radiation pressure |2|. 

We begin our analysis of equilibria of the electrodynamic Casimir force with 



the precursor of Eq. ( 39 1, which contains the abstract T and G^-operators, where 
Gm is the electromagnetic Green's function operator for an isotropic, homogeneous 
medium, [J 

Tic f 00 

S=— dK\x lnTr'TLo, (43) 

2n Jo 

where the operator [T _1 (icK,x,x')] equals 



4 Gm satisfies (V x n M l (icK)V x +£ M {icK)K 1 ) Gm(<ck,x,x') = 5(x - x')I, and is related to 
Gm, the Green's function of the imaginary frequency Helmholtz equation, by Gm('ck,x,x') = 
jIIm(z'cK') (l+ (nMK)~ 2 V(g> V') GM(icriMK,x,x'). Here, tiM(icK) = \J £M(ick)yiM(icK) is the index 
of refraction of the medium, whose argument is suppressed to simplify the presentation. Thus Gm, 
in contrast to Go, takes into account the permittivity and permeability of the medium when they 
are different from one. 
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[Tj 1 (icK,xi,x' 1 )] [G M (icK,x h x' 2 )} -\ 

[GM(icK,x 2 ,x\)} [Tg 1 (JcJc,x 2 ,4)] , (44) 

and Too is the inverse of T -1 with Gm set to zero. The square brackets "[ ]" denote 
the entire matrix or submatrix with rows indicated by x and columns by x'.jjThe 
operator [T _1 (icK,x,x')} has indices in position space. Each spatial index is limited 
to lie inside the objects A,B, ■ ■ ■ . For both indices x and x' in the same object A the 
operator is just the inverse T operator of that object, [T7 (/cfc,x,x')]. For indices on 
different objects, x in A and x' in B, it equals the electromagnetic Green's function 
operator [Gj}(i'cK 1 x,x')] for an isotropic, homogeneous medium. 

As shown in section |2.4| after a few manipulations, the operators Tj and Gm 
turn into the on-shell scattering amplitude matrix, ¥j, of object / and the translation 
matrix X, which converts wave functions between the origins of different objects. 
While practical computations require evaluation of the matrices in a particular wave 
function basis, the position space operators Tj and Gm are better suited to our gen- 
eral discussion here. 

To investigate the stability of object A, we group the 'rest' of the objects into a 



single entity R. So, T consists of 2 x 2 blocks, and the integrand in Eq. (43 1 reduces 
to tr In (I — TaGmTrGm)- Merging the components of R poses no conceptual dif- 
ficulty given that the operators are expressed in a position basis, while an actual 
computation of the force between A and R would remain a daunting task. If object 
A is moved infinitesimally by vector d, the Laplacian of the energy is given by 

+ 2T A VG M T R (VG M ) r i-T A GkG M ( 46 ) 

+ 2T A V<S M T R G M i _ JaGmTrGm (47) 

T A VG M T R G M + T A G M T R (WG M ) J j-t aGm j rGm 

After displacement of object A, the Green's function multiplied by T A on the left 
and ¥# on the right (T^GmTr) becomes GM(icK,X + d,x'), while that multiplied 
by Tr on the left and T A on the right (TrGmTa) becomes Gm(('c(C,x,x' -f-d). The 
two are related by transposition, and indicated by VGm(*'cK",x,x / ) = VdGjvf(jcJC,x+ 
d,x')|d = o and (VGm('ck - ,x,x / )) = VdGAf(/cK",x,x'-r-d)|d=o in the above equation. 
In the first line we have substituted n M {icK) k 2 Gm for V 2 Gm', the two differ only by 
derivatives of 5-functions which vanish since Gm (/cfc,x,x') is evaluated with x in 
one object and x' in another. In expressions not containing inverses of T-operators, 
we can extend the domain of all operators to the entire space: Tj(icK,x,x ! ) — if x 
or x' are not on object J and thus operator multiplication is unchanged. 



5 To obtain the free energy at finite temperature, in place of the ground state energy S, J || is 
replaced by the sum W- Yl K >o over Matsubara 'wavenumbers' K„ = lltnkT /he with the K"o = 
mode weighted by 1 /2. 
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To determine the signs of the various terms in Vjj<?|._ , an analysis similar to 
Ref. [51 1 can be performed. Consequently, the Laplacian of the energy is found 
to be smaller than or equal to zero as long as both T^ and Tg are either positive or 
negative semidefinite for all imaginary frequenciesH The eigenvalues of Tj, defined 



in Eq. ( 25 i, on the other hand, are greater or smaller than zero depending on the sign 
s J of Vj, since 

Tj = s j ^/?Vj 7 = ^=V^ 7 V/. (48) 

We are left to find the sign of the potential, 

Yj(icK,x) = Ik: 2 (ej(icK,x) - EmO'ck)) 

+ V x (fiy l (icK,x)-ix M l (icK)) Vx , (49) 

of the object A, and the compound object RVj The sign is determined by the rela- 
tive permittivities and permeabilities of the objects and the medium: If Ej(icK,x) > 
£m{ick) and jj.j(icK,x) < plm{ick) hold for all x in object 7, the potential Vj is posi- 
tive. If the opposite inequalities are true, Yj is negative. The curl operators surround- 
ing the magnetic permeability do not influence the sign, as in computing an inner 
product with Vy they act symmetrically on both sides. For vacuum Em = Mm = 1> 
and material response functions e(icK,x) and jj,(icK,x) are analytical continuations 
of the permittivity and permeability for real frequencies ll60ll . While E(icK,x) > 1 
for positive K, there are no restrictions other than positivity on jj,(icK,x). (For non- 
local and non-isotropic response, various inequalities must be generalized to the 
tensorial operators £ (icK",x,x') and jtl (icK,x,x').) 

Thus, levitation is not possible for collections of objects characterized by Ej(icK,x) 
and jj.j(icK,x) falling into one of the two classes described earlier, i) £j/Em > 1 and 
1-i-j/fJ-M < 1 (positive Yj and Tj), or ii) Ej/Em < 1 and \IjJ\Im > 1 with (negative 
Vj and Tj). (Under these conditions parallel slabs attract.) The frequency and space 
dependence of the functions has been suppressed in these inequalities. In vacuum, 
£m{ick) — jj,M{icK) — 1; since E(icK,x) > 1 and the magnetic response of ordi- 
nary materials is typically negligible |60|, one concludes that stable equilibria of 
the Casimir force do not exist. If objects A and R, however, belong to different cat- 
egories — under which conditions the parallel plate force is repulsive — then the 
terms under the trace in lines ([45} and (|46| are negative. The positive term in line 



(47 1 is typically smaller than the first two, as it involves higher powers of T and 
Gm- In this case stable equilibrium is possible, as demonstrated recently for a small 
inclusion within a dielectric filled cavity [79]. For the remaining two combinations 
of inequalities involving £jJ£m and flj/flM the sign of Yj cannot be determined a 
priori. But for realistic distances between objects and the corresponding frequency 



6 In practice, Tx and T^ suffice to have the same sign over the frequencies, which contribute most 
to the integral (or the sum) in Eq. d43b . 

7 The first curl in the operator Yj results from an integration by parts. It is understood that it acts 
on the wave function multiplying V/ from the left. 
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ranges, the magnetic susceptibility is negligible for ordinary materials, and the in- 
equalities involving ji can be ignored. 

In summary, the instability theorem applies to all cases where the coupling of 
the EM field to matter can be described by response functions e and jj., which may 
vary continuously with position and frequency. Obviously, for materials which at a 
microscopic level cannot be described by such response functions, e.g., because of 
magneto-electric coupling, our theorem is not applicable. 

Even complicated arrangements of materials obeying the above conditions are 
subject to the instability constraint. For example, metamaterials, incorporating ar- 
rays of micro-engineered circuity mimic, at certain frequencies, a strong magnetic 
response, and have been discussed as candidates for Casimir repulsion across vac- 
uum. (References [87, 86 1 critique repulsion from dielectric/metallic based meta- 
materials, in line with our following arguments.) In our treatment, in accord with 
the usual electrodynamics of macroscopic media, the materials are characterized 
by e(icK,x) and fi(icK,x) at mesoscopic scales. In particular, chirality and large 
magnetic response in metamaterials are achieved by patterns made from ordinary 
metals and dielectrics with well-behaved e(icK,x) and fi(icK,x) sa 1 at short scales. 
The interesting EM responses merely appear when viewed as 'effective' or 'coarse 
grained'. 

Clearly, the coarse-grained response functions, which are conventionally em- 
ployed to describe metamaterials, should produce, in their region of validity, the 
same scattering amplitudes as the detailed mesoscopic description. Consequently, 
as long as the metamaterial can be described by e(icK,x) and jj,(icK,x) w 1, the 
eigenvalues of the T operators are constrained as described above, and hence sub- 
ject to the instability theorem. Thus, the proposed use of chiral metamaterials in 
reference [96] cannot lead to stable equilibrium since the structures are composites 
of metals and dielectrics. Finally, we note that instability also excludes repulsion 
between two objects that obey the above conditions, if one of them is an infinite 
flat plate with continuous translational symmetry: Repulsion would require that the 
energy as a function of separation from the slab should have d^S > at some point 
since the force has to vanish at infinite separation. A metamaterial does not have 
continuous translational symmetry at short length scales but this symmetry is ap- 
proximately valid in the limit of large separations (long wavelengths), where the 
material can be effectively described as a homogeneous medium. At short separa- 
tions lateral displacements might lead to repulsion that, however, must be compati- 
ble with the absence of stable equilibirum. 



4 Applications 

This section gives an overview on different geometries and shapes that have been 
studied by the approach that we introduced in Section[2] A selection of applications 
has been made to showcase generic situations and important effects that had not 
been studied in detail before the development of the methods described here. We 
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shall mainly summarize analytical and numerical results for the Casimir interaction 
in the various systems. For details on their derivation and additional implementa- 
tions of the scattering approach we refer to the literature. 



4.1 Cylinders, wires, and plates 

The extent to which EM field fluctuations are correlated depends on the effective 
dimensionality of the space that can be explored by the fluctuations. Therefore, 
Casimir interactions are expected to depend strongly on the codimension of the 
interacting objects. The focus of this subsection is on the particular properties of 
systems with a codimension of the critical value two. We consider these problems 
in the context of interactions between cylinders and a cylinder and a plate, both 
perfect reflectors and dielectric materials. Cylindrical geometries are of recent ex- 
perimental interest since they are easier to hold parallel than plates and still generate 
a force that is extensive in one direction. 



Fig. 5 (a) Casimir energy for 
two cylinders of equal radius 
R as a function of surface- 
to-surface distance d — 2R 
(normalized by the radius). 
The energy is divided by the 
PFA estimate E^p^ for 
the energy. The solid curves 
show our numerical results; 
the dashed lines represent 
the asymptotic results of 
Eq. {54}. (b) Casimir energy 
for a cylinder of radius R 
parallel to a plate as a function 
of the surface-to-surface 
distance H — R (normalized 
by the radius). The energy is 
divided by the PFA estimate 
E L p FA p e . The solid curves 
reflect our numerical results; 
the dashed lines represent the 
asymptotic results of Eq. l|58[l. 



(a) 
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We consider two cylinders of equal radii R and length L — > °° with center-to- 
center separation d, see Fig.[5ta) [76 1. (The related configuration where one cylin- 
der is inside another cylinder is treated in reference [22 j.) For this geometry the 
interaction energy is obtained from the expression 

g= T r dKl °% det (f - ^!y\^ ba ^!y\^ ab ) ■ (50) 



with the exterior scattering amplitudes of a cylinder. 

t ^cyl,k z n'E,k z nM — ^cyl,k z n'M,k z nE " 



^cyl,k z ii'M ; k z nM ~ L°V- z ' ' K' (Rd) ' (51) 



^cyl,k z n'E,k z nE — l"^ z ' ' K (Rp) 



and the matrices a tt ab , °?/ ba that translate from cylinder a to b and vice versa. Their 
elements are summarized in Ref. [75]. The matrix inside the determinant is diag- 
onal in k z , so the log-determinant over this index turns into an overall integral. A 
change of variable to polar coordinates converts the integrals over K and k z to a 
single integral over p = yjk 2 + K 2 , yielding 

/ pt/p(logdet^' M + logdet^ /£ ), (52) 



where 



47T Jo 



^n.n" = 5 "-""-\Y^pM) Kn+Apd) ^A^) n ' +n " {P) 



(53) 



describe magnetic (TE) or Neumann modes and electric (TM) or Dirichlet modes, 
respectively. 

For large separations d^> R, the asymptotic behavior of the energy is determined 
by the matrix elements for n = nf = for Dirichlet modes and n = n' = 0, ±1 for 
Neumann modes. Taking the determinant of the matrix that consists only of these 
matrix elements and integrating over p yields straightforwardly the attractive inter- 
action energies 

**=-?£ ! ii- 



d 2 8nlog 2 (d/R) \ \og{d/R) ■ ,. (>1) 

gM^_ hcL l_RL_ 
5nd 6 
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for electric (Dirichlet) and magnetic (Neumann) modes. The asymptotic interaction 
is dominated by the contribution from electric (Dirichlet) modes that vanishes for 
R — >• only logarithmically. 

For arbitrary separations higher order partial waves have to be considered. The 
number of partial waves has to be increased with decreasing separation. A numeri- 
cal evaluation of the determinant and the /^-integration can be performed easily and 
reveals an exponentially fast convergence of the energy in the truncation order for 
the partial waves. Down to small surface-to-surface separations of (d — 2R)/R = 0.1 
we find that n = 40 partial waves are sufficient to obtain precise results for the en- 
ergy. The corresponding result for the energies of two cylinders of equal radius is 
shown in Fig. |5ja). Notice that the minimum in the curve for the total electromag- 
netic energy results from the scaling by the proximity force approximation (PFA) 
estimate of the energy. The total energy is monotonic and the force attractive at all 
separations. 

Next we consider a cylinder and an infinite plate, both perfectly reflecting, see 
Fig. |5jb). The Casimir energy for this geometry has been computed originally in 
Ref. [33|. In the limit of perfectly reflecting surfaces, the method of images can be 
employed to compute the Casimir interaction for this geometry [76 1. Here we use a 
different method that can be also applied to real metals or general dielectrics ll75ll . 
We express the scattering amplitude of the cylinder now in a plane wave basis, using 

^cyl.kj_P.k' P' _ 2- r U k ± P,k,nQ^cyLk z nQ,kzn'Q' U k z n>Q>.k' ± P'> ^J 

nQ,n'Q' L 2 

where kx denotes the vector (k y ,k z ), C^, p(k) and Cq are normalization coefficients 
that can be found together with the matrix elements of the conversion matrix D in 
Ref. 11751 . The elements of the scattering amplitude in the cylindrical basis are given 
by Eq. (J5T|i. The scattering amplitude of the plate is easily expressed in the plane 

wave basis as 

^plate,k' E,k ± M — "^ plate X ' M,kj_E ~ u ' 



^ate,k' ± M.k ± M = ^^ (k ± - k ± ) r» [icK, ^1 + tfjtf j. (%) 

^o7ate.k',E.k,E = ^^(kx-k'Jr* ( icK, J I + k 2 jK 2 ' 



plate, k^E.k^E 

in terms of the Fresnel coefficients that read for a general dielectric surface 

^ = n{icK)-^l + (n 2 (icK)-l)^ 

|l(icif) + ;i|(fl 2 (M)-l)i 2 ' 

p.. . e(icK) -J\ + (n 2 (icK) - \)x 2 

r-(icK,x) = — v . 

e(icK) + VI + (n 2 (icK) - \)x 2 
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Here, n is the index of refraction, n{icK) = ^e{icK)n{icK). In the limit of a per- 
fectly reflecting plate one has r M — >• — 1, r E — >• 1. The energy given by Eq. (50 1 can 



now be evaluated in the plane wave basis with the translation matrices given by the 

simple expression ^ pk , p , = e _ V k ^ ff ^ 2 5( 2 )(k_ L -* x )8pp. 

The asymptotic expression for the attractive interaction energy at large distance 
// > R reads 

„ F hcL 1 



H* 167t]Qg(H/R)' (5g) 

32k H 4 

The total electromagnetic Casimir interaction is again dominated by the contribution 
from the electric (Dirichlet) mode with n = which depends only logarithmically 
on the cylinder radius. The interaction at all separations follows, as in the case of 



two cylinders, from a numerical computation of the determinant of Eq. ( 50 » and 
integration over p. The result is shown in Fig.[5jb). 

The above approach has the advantage that it can be also applied to dielectric 
objects. The scattering amplitude of a dielectric cylinder can be obtained by solving 
the wave equation in a cylindrical basis with appropriate continuity conditions ||75ll . 
The scattering amplitude is diagonal in k z and the cylindrical wave index n, but not 
in the polarization. Here we focus on large distances H S> R- Expanding the logdet 



in Eq. (50 1, we obtain for the interaction energy 

S = - lUS t dX ^V\ [(? + £cyl '° " 4 *> £ (°^ - ( 3 + £cyi,o)* 2 r>,*)] , 

(59) 
if the zero-frequency magnetic permeability jJ. cy ifi of the cylinder is set to one. If we 
do not set /i cy i,o equal to one, but instead take the perfect reflectivity limit for the 
plate, we obtain 

_ _ hcLR 2 (£ C yl,0 ~ Mcyl.o) (9 + £cyl,0 + Mcyl.O + £cyl,OMcyl,o) ™ 

327T// 4 ~ (l+£cyl,0)(l+Mcyl,0) "' 

Finally, if we let £ cy i be infinite from the beginning (the perfect metal limit for 
the cylinder), only the n = TM mode of the scattering amplitude contributes at 
lowest order. For a plate with zero-frequency permittivity £ p i a te,o and permeability 
Mpiate,o> we obtain for the Casimir energy 

S= l6nJ\o g {R/H/^ m 

where 

-1 dx 



t E = / ^Mr^O,*)-*^,*)]. (62) 

Jo l+x L 
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In Fig.pl £ is plotted as a function of the zero-frequency permittivity of the plate, 
£piate.Cb f or various zero-frequency permeability values, jU p i a te,o- 



Fig. 6 Plots of <p E versus 
l/fiplate.o for fixed values of 
A^plate.o- The perfect metal 
limit (0 £ = 1) is approached 
slowly for large /i p iate,o> a s in 
the case of a sphere opposite 
a plate. For large /ipi ate ,o the 
interaction becomes repul- 
sive, which is expected given 
similar results for two infinite 
plates. 




4.2 Three-body effects 



Casimir interactions are not pair-wise additive. To study the consequences of this 
property, we consider the case of two identical objects near perfectly reflecting 
walls Il78ll76ll . Multibody effects were first observed for such a configuration with 
two rectangular cylinders sandwiched between two infinite plates by Rodriguez et 
al. 1 83 1. The role of dimension on this effect is studied by considering either cylin- 
ders, see Fig. 17] or spheres, see Fig. [8] While we have given a more detailed de- 
scription of how the interaction energies follow from the scattering approach in the 
previous subsection, we mainly provide the final results in this and in the following 
subsections. 

First, we consider the geometry shown in Fig.|7jwith two cylinders that are placed 
parallel to one or in-between two parallel plates, where all objects are assumed to 
be perfectly reflecting. Using the general expression for the Casimir energy of mul- 
tiple objects, Eq. ( f39] l, the energy can be straightforwardly computed by truncating 
the matrix M at a finite partial wave order n. Including up to n = 35 partial waves, 
we obtain for the Casimir force between two cylinders of equal radii in the pres- 
ence of one or two sidewalls the results shown in Fig. U\ In this figure the force at 
a fixed surface-to-surface distance d — 2R = 2R between the cylinders is plotted as 
a function of the relative separation (H — R)/R between the plate and cylinder sur- 
faces. Two interesting features can be observed. First, the attractive total force varies 
non-monotonically with H: Decreasing for small H and then increasing towards the 
asymptotic limit between two isolated cylinders for large H, cf. Eq. ( |54| . The ex- 
ternum for the one-sidewall case occurs at H — R ss 0.27 R, and for the two-sidewall 
case is at H — R « QA6R. Second, the total force for the two-sidewall case in the 
proximity limit H = R is larger than for H/R — > °°. As might be expected, the H- 
dependence for one sidewall is weaker than for two sidewalls, and the effects of the 
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Fig. 7 Electromagnetic Casimir force per unit length between two cylinders of radius R and side- 
wall separation H vs. the ratio of sidewall separation to cylinder radius (H — R)/R, at fixed distance 
(d — 2R)/R = 2 between the cylinders, normalized by the total PFA force per unit length between 
two isolated cylinders, Fpfa = | (/icTT 3 / 1920) \J R / (d — 2R) 1 . The force is attractive. The solid 
lines refer to the case with one sidewall, while dashed lines depict the results for two sidewalls. 
Also shown are the individual TE (blue) and TM (red) forces. 



two sidewalls are not additive: not only is the difference from the H — > °° force not 
doubled for two sidewalls compared to one, but the two curves actually intersect at 
a separation of H/R = 1.13. The non-monotonic sidewall effect arises from a com- 
petition between the force from TE and TM modes as demonstrated by the results in 
Fig. [7] The qualitatively different behavior of TE and TM modes can be understood 
intuitively on the basis of the method of images [ 76 1. The non-monotonicity in H 
also implies that the force between the cylinders and the sidewalls is not monotonic 
in d [76). 

Second, we replace the two cylinders by two identical, general polarizable com- 
pact objects that we specialize later on to spheres |85|. The meaning of the lengths 
d and H remains unchanged. In dipole approximation, the retarded limit of the in- 
teraction is described by the static electric (a z , a«) and magnetic (/3,, j3ii) dipole 
polarizabilities of the objects which can be different in the directions perpendicu- 
lar (z) and parallel (||) to the wall. The well-known Casimir-Polder (CP) potential 
between two compact objects at large distance is 



\(d) = - 



lie 
87td 7 



33a|? + 13a 7 2 -14aiij3 z + (aoj3) 



(63) 



When a sidewall is adde , the energy changes. Its ^/-dependent part is then 



<&o (d,H) = # 2> | id) + <? 2 ,\ (D,d) + S-i (D, d) 



(64) 
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22 : F/F(H -> oo 



1.0 



Dots top to bottom: 

R/d = 0.05, 0.1 0,0.1 5, 

0.20,0.25,0.30 




1 2 (H - R)/R 3 4 5 

Fig. 8 Electromagnetic Casimir force between two spheres next to one sidewall at separation H 
vs. the ratio H/R for different sphere separations d. Dotted curves represent numerical results. 
Shown are also the analytical results of Eq. (68) , including terms up to j = 10 for R/d < 0.2 (solid 
curves) |85|. Inset: Magnification of the nonmonotonicity. 



with D = \Jd 2 +4H 2 . The change in the relative orientation of the objects with 
£ = d/D leads to a modification of the 2-body CP potential 



&\(D,d) = 



he 



87TD 7 

63£ 4 (a,i 



■a z ) -14(a,|/3||(l 



2 "of 

2\ , ,J 



26a,f +20o£ - 14T (4a,f - 9ai a, + 5a£) 



«llj3 z ) + (a^/3)] 



(65) 



The three-body energy $$ (D, d) describes the collective interaction between the two 
objects and one image object. It is given by 



i(p,d) 



Ahc 



1 



% d 3 D 4 {£+l) 5 
xfaj-pfj-fol 6 

+ 4(£ 4 + 5£ 3 +f) (a z P\\-a\\P z 



3F + 15£ 3 + 2%t + 20£ 3 + 6£ l - 5£ - 1 
15^ 5 +24^ 4 - \Q£ 2 - 5£ - l) (a? - j8*) 



(66) 



It is instructive to consider the two limits H <C d and H ^ d. For H <</, <? 00 turns 
out to be the CP potential of Eq. ( |63] l with the replacements a, — » 2a z , (Xii — > 0, 
j3 z — > 0, j8ii — » 2j8ii . The two-body and three-body contributions add constructively or 
destructively, depending on the relative orientation of a dipole and its image which 
together form a dipole of zero or twice the original strength l85l . 
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For H~^> d the leading correction to the CP potential of Eq. ( |63| i comes from the 
three-body energy. The energy then becomes (up to order H~ 6 ) 



he 
g ao {d,H) = £ 2 \{d) + — 

— " n 



a? — ajf 9a,? - a? - 2a\S z 



(67) 



The signs of the polarizabilities in the leading term ~ H~ 4 can be understood from 
the relative orientation of the dipole of one object and the image dipole of the other 
object |85|. 

Next, we study the case where the two objects are perfectly reflecting spheres of 
radius R. Now we consider arbitrary distances and include higher order multipole 
contributions. For R-^id,H and arbitrary H/d the result for the force can be written 
as 

The functions /,• can be computed exactly and their full form is given for j — 
6,7, 8 in Ref. [85|. For H > d one has f 6 (h) = -1001/16 + 3/(4/i 6 ) + ^(/r 8 ), 
h(h) = -71523/160 + 39/(80/i 6 ) + ^(/T 8 ) so that the wall induces weak re- 
pulsive corrections. For H < d, f 6 (h) = -791/8 + 674lh 2 /8 + 0(h 4 ), fy{h) = 
-60939/80 + 582879/i 2 /80 + 0(h A ) so that the force amplitude decreases when 
the spheres are moved a small distance away from the wall. This proves the exis- 
tence of a minimum in the force amplitude as a function of H for fixed, sufficiently 
small R/d. 

To obtain the interaction at smaller separations or larger radius, the energy S aa 
and force F — —dS'oo/dd between the spheres has been computed numerically l85ll . 
In order to show the effect of the side wall, the energy and force between the spheres, 
normalized to the results for two spheres without a wall, is shown in Fig. [8] for 
fixed d. When the spheres approach the wall, the force first decreases slightly if 
R/d < 0.3 and then increases strongly under a further reduction of//. For R/d > 0.3 
the force increases monotonically as the spheres approach the wall. This agrees with 
the prediction of the large distance expansion. The expansion of Eq. ( |6"8"j ) with j = 10 
terms is also shown in Fig.jslfor R/d < 0.2. Its validity is limited to large d/R and 
not too small H/R; it fails completely for R/d > 0.2 and hence is not shown in this 
range. 



4.3 Orientation dependence 

In this subsection we describe the shape and orientation dependence of the Casimir 
force using Eq. ([39|, first reported in Ref. 1 30]. We consider the orientation depen- 
dent force between two spheroids, and between a spheroid and a plane. For two 
anisotropic objects, the CP potential of Eq. ( |63"] > must be generalized. In terms of the 
Cartesian components of the standard electric (magnetic) polarizability matrix a 
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(j3), the asymptotic large distance potential of two objects (with the z axis pointing 
from one object to the other), can be written as 



S = 



- 20 alal z - 30 (c^a 2 , + 0,1a 2 ,) + (a ->• j8 ) 
7 (a.4/3 2 , + ^2 _ 2^2 ) + (1 ^ 2) | . 



(69) 



For the case of an ellipsoidal object with static electric permittivity e and magnetic 
permeability ju, the polarizability tensors are diagonal in a basis oriented to its prin- 
cipal axes, with elements (for i € {1,2,3}) 



"' An 1 + (e - 1 )», ' "* An 1 + Qu - l)n t 



(70) 



where V = 4nr\r2ri/3 is the ellipsoid's volume. In the case of spheroids, for which 
r\= r2=R and rj, = L/2, the so-called depolarizing factors, n ; -, can be expressed in 

2e), where the 



terms of elementary functions, n\ = «2 = ^r 1 : n 3 = ^^(log fz 



l+e 
1- 



eccentricity e = (1 — %*- ) ' is real for a prolate spheroid (L > 2/?) and imaginary for 
an oblate spheroid (L < 2R). The polarizability tensors for an arbitrary orientation 
are then obtained as a — S%~ l oPS%, where Stf, is the matrix that orients the prin- 
cipal axis of the spheroid relative to a fixed Cartesian basis. Note that for rarefied 
media with £ ~ 1, fl ~ 1 the polarizabilities are isotropic and proportional to the 
volume. Hence, to leading order in e — 1 the interaction is orientation independent 
at asymptotically large separations, as we would expect, since pairwise summation 
is valid for e — 1 <C 1 . In the following we focus on the interesting opposite limit of 
two identical perfectly reflecting spheroids. We first consider prolate spheroids with 
L^> R. The orientation of each "needle" relative to the line joining them (the initial 
z-axis) is parameterized by the two angles (9, \jf), as depicted in Fig.[9ta). Then the 
energy is 



he f 5L 6 

« 1 115271 (in |-l) 



1 7 

COS 01 COS 02 



+ — cos 2 v/ sin 2 0i sin 2 02 - -cos y/ sin 2 0i sin 2 02 
20 8 



L*R l 



(71) 



where \]/= y/i — y/2- It is minimized for two needles aligned parallel to their separa- 
tion vector. At almost all orientations the energy scales as L 6 , and vanishes logarith- 
mically slowly as R — > 0. The latter scaling changes when one needle is orthogonal 
to z (i.e. 0i = n/2), while the other is either parallel to z (02 = 0) or has an arbitrary 
02 but differs by an angle n/2 in its rotation about the z-axis (i.e. y/\ — 1//2 = n/2). 



In these cases the energy comes from the next order term in Eq. (71 1, and takes the 
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form 

/ (f-«-f)- i^4-» ca+7c " 2 * ) - <72) 

K 

which shows that the least favorable configuration corresponds to two needles or- 
thogonal to each other and to the line joining them. 

For perfectly reflecting oblate spheroids with R ^> L/2, the orientation of each 
"pancake" is again described by a pair of angles (9, y), as depicted in Fig.|9jb). To 
leading order at large separations, the energy is given by 



S- 



he f R 6 

^\l447T 3 



765-5(cos20 1 +cos20 2 )+237cos20 1 cos20 2 

$(R 5 L) 



- 372cos2y/sin 2 9\ sin 2 2 — 300cos y/sin20i sin202 



(73) 



The leading dependence is proportional to R 6 , and does not disappear for any choice 
of orientations. Furthermore, this dependence remains even as the thickness of the 
pancake is taken to zero (L — > 0). This is very different from the case of the needles, 
where the interaction energy vanishes with thickness as ln~ l (L/R). The lack of L 
dependence is due to the assumed perfectly reflectivity. The energy is minimal for 
two pancakes lying on the same plane (0i = 62 = 7C/2, y/ = 0) and has energy 
— he {\12> / \%n i )R 6 / d 1 . When the two pancakes are stacked on top of each other, 
the energy is increased to —hc{62/9K i )R 6 /d 1 . The least favorable configuration is 
when the pancakes lie in perpendicular planes, i.e., Q\ = n/2, O2 = 0, with an energy 
-he {11 /3n 3 )R 6 /d 1 . 

For an anisotropic object interacting with a perfectly reflecting mirror, at leading 
order the CP potential is given by 

Pie 1 

S = -- A -tr{a-^) + e{d- 5 ) 1 (74) 

which is clearly independent of orientation. Orientation dependence in this system 
thus comes from higher multipoles. The next order also vanishes, so the leading 
term is the contribution from the partial waves with I = 3 for which the scattering 
matrix is not known analytically. However, we can obtain the preferred orientation 
by considering a distorted sphere in which the radius R is deformed to R + 5/(#, <p) . 
The function / can be expanded into spherical harmonics Y( m (d-,<p), and spheroidal 
symmetry can be mimicked by choosing / = Y2o(#,<p)- The leading orientation 
dependent part of the energy is then obtained as 

S f = -hc 16 y * cos(28). (75) 

A prolate spheroid (8 > 0) thus minimizes its energy by pointing towards the mirror, 
while an oblate spheroid (8 < 0) prefers to lie in a plane perpendicular to the mirror. 
(It is assumed that the perturbative results are not changed for large distortions.) 
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Fig. 9 (a) Orientation of 
a prolate (cigar-shaped) 
spheroid: The symmetry axis 
(initially the z-axis) is rotated 
by 8 about the x-axis and then 
by \j/ about the z-axis. For two 
such spheroids, the energy 
at large distances is give by 
Eq. l |7H . The latter is depicted 
at fixed distance d, and for 
Vi — ¥2 , by a contour plot as 
function of the angles 0| , 02 
for the x-axis rotations . Min- 
ima (maxima) are marked by 
filled (open) dots, (b) As in (a) 
for oblate (pancake-shaped) 
spheroids, with a contour plot 
of energy at large separations. 




These configurations are also preferred at small distances d, since (at fixed distance 
to the center) the object reorients to minimize the closest separation. Interestingly, 
the latter conclusion is not generally true. In Ref. If30ll it has been shown that there 
can be a transition in preferred orientation as a function of d in the simpler case 
of a scalar field with Neumann boundary conditions. The separation at which this 
transition occurs varies with the spheroid's eccentricity. 



4.4 Edge and finite size effects 



In this subsection, based on work reported in Ref. [45|, it is demonstrated that 
parabolic cylinders provide another example were the scattering amplitudes can 
be computed exactly. We use the exact results for scattering from perfect mirrors to 
compute the Casimir force between a parabolic cylinder and a plate. In the limiting 
case when the curvature at its tip vanishes, the parabolic cylinder becomes a semi- 
infinite plate (a knife's edge), and we can consider how edges and finite size effects 
influence the Casimir energy. 

The surface of a parabolic cylinder in Cartesian coordinates is described by y = 
(x 2 — R 2 )/2R for all z, as shown in Fig. 10 a), where R is the curvature at the tip. In 
parabolic cylinder coordinates, defined through x — JJ.X, y — (X 2 — /i 2 )/2, z = z, the 
surface is simply jj,= jUo — >/R for — °o < X , z < °°. Since sending X — >• — X and /I — > 
— ju returns us to the same point, we restrict our attention to fi > while considering 
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all values of A. Then fi plays the role of the "radial" coordinate in scattering theory 
and one can again define regular and outgoing waves P31 . Since both objects are 
perfect mirrors, translational symmetry along the z-axis enables us to decompose 
the electromagnetic field into two scalar fields, as in the case of circular cylinders in 



Subsection 4.1 Each scalar field, describing E (Dirichlet) orM (Neumann) modes, 
can then be treated independently, with the sum of their contributions giving the full 
electromagnetic result. 

The scattering amplitude of the plate is expressed in a plane wave basis and 
is given by Eq. ( |56| > with r 14 = — 1 and r E = 1. The scattering amplitude of the 
parabolic cylinder for E and M polarization is obtained in a parabolic cylinder wave 
basis as 1451 



(76) 



aree 27t xn tJ\X t r _ -v D v(ifk) 

^para,k z vE,k£v'E-~^ (A-/y<) v y/ <:cV £, J k _ vE - I — — — 

,&ee _ 2n x(, ,t\x f f -V+l D v( l Po) 

^para,k z vM,k z v'M ~ ~ -^ { K Z ~ K z) vyJk z vM, Jk z vM - ' Jy 7~-y , 



with fiQ = JlR^k 1 +k} and the parabolic cylinder function D v (u) for integer v. 
For the present geometry, the general formula for the Casimir energy per unit 
length can be expressed explicitly as 

^l = Io^I-^ log det ( dvy ~ fk - yp I-i kx ® vkxkz ( - d ' e)rF ^ v ' k ^ {d ~ e 

(77) 
for polarization P = E or M. Here the matrix ^ with elements 



— (tan^r 
lkyV\y/2% cos- 



%k xk -M,0) = X L,_..: ^ X -ij-^ W 



with k y = i^/k 2 +k^+k^ and tan = k x /k y describes the translation from parabolic 
cylinder to plane waves over the distance d from the focus of the parabola to the 
plane where 9 is the angle of inclination of the parabolic cylinder. 

Numerical computations of the energy are performed by truncating the determi- 
nant at index v max . For the numbers quoted below, we have computed for v ram up to 



200 and then extrapolated the result for v milI — > °°, and in Fig. 10 we have generally 



used v mux = 100. The dependence of the energy on the separation H = d — R/2 for 



8 = is shown in Fig. 10 a). At small separations (H/R <C 1) the proximity force 
approximation, given by 

%1 r dX ■ (79) 



HcL 120j-~[H+x 2 /(2R)) 3 960n/2VH 5 ' 
should be valid. The numerical results in Fig. fTOja) indeed confirm this expectation. 
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Fig. 10 (a) Energy 
EH 2 /{hcL) versus H/R for 
= and R = 1 on a log- 
linear scale for the parabolic 
cylinder-plane geometry. The 
dashed line gives the R = 
limit and the solid curve gives 
the PFA result, (b) The coef- 
ficient c(6) as a function of 
angle for R = 0. The exact 
result at 9 = it/1 is marked 
with a cross. Inset: Dirich- 
let (circles) and Neumann 
(squares) contributions to the 
full electromagnetic result. 
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A more interesting limit is obtained when R/H —} 0, corresponding to a semi- 
infinite plate for which the PFA energy vanishes. The exact result for the energy for 

R = and = is 

S _ _ C± 



(80) 



where C± — 0.0067415 is obtained by numerical integration. When the semi-infinite 
plate is tilted by an angle 0, dimensional analysis suggest for the Casimir energy 



ficL 



C(B) 
H 2 



(81) 



The function c(9) = cos(9)C(0) is shown in Fig. 10 b). A particularly interesting 



limit is 9 — >• k/2, when the two plates are parallel. In this case, the leading con- 
tribution to the Casimir energy should be proportional to the area of the half -plane 
according to the parallel plate formula, Eu/{hcA) = — en/// 3 with cii = K 2 /720, 
plus a subleading correction due to the edge. Multiplying by cos 6 removes the di- 
vergence in the amplitude C(6) as 6 — > k/2. As in |39|, we assume c(6 — > k/2) = 
cii /2 + (0 — k/2) c edge , although we cannot rule out the possibility of additional non- 
analytic forms, such as logarithmic or other singularities. With this assumption, we 
can estimate the edge correction c LllgL . = 0.0009 from the data in Fig. 10 b). From 
the inset in Fig.[T0|b), we estimate the Dirichlet and Neumann contributions to this 
result to be c~L = —0.0025 and c^ c = 0.0034, respectively. For extensions to other 



Geometry and material effects in Casimir physics - Scattering theory 37 

geometries with edges, inclusion of thermal fluctuations and experimental implica- 
tions, see Ref. |45l . 



4.5 Interior configurations 

In this last subsection we consider so-called interior configurations where one ob- 
ject is contained within another that can be also studied with the methods introduced 
in Sect. [2] Specifically, we obtain the electrodynamic Casimir interaction of a con- 
ducting or dielectric object inside a perfectly conducting spherical cavity [95]. In the 
case where an object, i, lies inside a perfectly conducting cavity, the outer object o, 



the Casimir energy of Eq. ( 42 1 becomes 



2% h s det(^-J%'J^ e ) ' v ' 

where J^" is the scattering amplitude for interior scattering of the conducting cav- 
ity, a sphere in our case, and ^f e the scattering amplitude of the interior object. 
The amplitude matrix for interior scattering is the inverse of the corresponding ex- 
terior matrix. These scattering amplitudes are evaluated in a spherical vector wave 
basis with respect to appropriately chosen origins within each object. The transla- 
tion matrices, W l ° and r f w , relate regular wave functions between the coordinate 
systems of the interior object and the spherical cavity, see Ref. [75] for details. The 



determinant in the denominator of Eq. ( 82 1 subtracts the Casimir energy when the 



origins of the two objects coincide. This way of normalizing the Casimir energy 
differs from the exterior cases considered before, where the objects are removed to 
infinite separation; a choice that would be unnatural in the interior case. 

First, we determine the forces and torques on a small object, dielectric or con- 
ducting, well separated from the cavity walls. This is the interior analogue of the 
famous Casimir-Polder force on a polarizable molecule near a perfectly conducting 
plate [15|. In this case the first term in a multiple scattering expansion, where the 



integrand of Eq. (82 1 is replaced by —Tx{,'^"W ^f e 'i / '"), already gives an excel- 
lent approximation to the energy. Since the object is small, the scattering amplitude 
^ilmPl'mfP" ( wnere I ar, d m m ' e angular momentum indices and P labels M or E po- 
larization) can be expanded in powers of K. Only the following terms contribute to 
lowest order: &f e . „, ,Jk) = 2jc 3 0f£ ,/3 + 0(k a ), where a p , is the static electric 

i,\mP,lm'P\ ' mm * ^ >' mm 

(P = E) or magnetic (P — M) polarizability tensor of the inner object. We consider 
an exterior spherical shell of radius R and define a to be the displacement of the 
center of the interior object from the center of the shell. Using the dipole approxi- 
mation for the inner object but including all multipoles of the exterior shell, we find 
for the Casimir energy to leading order in r/R (where r is the typical length scale of 
the interior object), the energy 
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(83) 



The z-axis is oriented from the center of the shell to the innterior object, and af: 
represent the interior object's static polarizability tensors in a Cartesian basis. The 
coefficient functions f p and g p can be obtained in terms of an integral over modified 
Bessel functions, see Ref. 11951 . f E is negative and decreasing with a/R, while f M 



is positive and increasing. There are important differences between Eq. (83 i and the 
classic Casimir-Polder result: first, the energy depends in a non-trivial way on a/R; 
second, at any non-zero distance from the center, the interior object experiences a 
torque; and third, the force between the two bodies depends on the interior object's 
orientation. 



To explore the orientation dependence of Eq. ( 83 i assume, for simplicity, there 
is a single frame in which both a E and a M are diagonal. In this body-fixed frame, 
write ctx X — a° v = j3 and aZ - \ (a® x + a&) = y (where we have suppressed the M/E 
label). The polarizability in the "lab frame" is obtained by a = &Op3%~ , where ££ 
is a rotation matrix that orients the principal axes of the inner object with respect to 
the lab frame. This procedure leaves Tra° invariant, and gives for the second line in 
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Eq. (83), 

£ /(a/fi) ( — sin 2 0cos2^ + 7 /> (3cos 2 0-i; 

P=M,E V 2 

where corresponds to the azimuthal rotation of the object about its principal z- 
axis, and 9 is the angle between the object's principal z-axis and the "laboratory" 
z-axis connecting the center of the sphere to the origin of the interior object. 

If /3 t^ then the object held at fixed inclination, 9, experiences a torque that 
causes it to rotate about the body-fixed z-axis. If, however, the object has axial sym- 
metry (j3 = 0), then the only torque on the object tries to align it either parallel or 
perpendicular to the displacement axis. 

A "cigar shaped" object (y > 0) prefers to orient so as to point perpendicular to 
the z axis, and a "pancake" (y < 0) tries to align its two large axes perpendicular to 
the z axis. The small ellipse inside the sphere in Fig. HIT a) illustrates a side view of 
both the cigar and the pancake in their preferred orientation. It is interesting to note 
that g E and g^ 1 are both positive. So, in contrast to the force, the contributions to the 
torque from magnetic and electric polarizabilities are in the same direction, if they 
have the same sign. More complicated behavior is possible if, for example, the elec- 
tric and magnetic polarizabilities are not diagonal in the same body-fixed coordinate 
system. Note that our results cannot be compared to the PFA approximation since 
the the size of the inner object, not the separation of surfaces, d, has been assumed 
to be the smallest scale in the analysis. 

An identical analysis can be performed for a polarizable object outside a metallic 
sphere where a/R > 1. It turns out that the analogous exterior function g(a/R) < 
for both polarizations. Therefore, the preferred orientation of a polarizable object 
outside a metallic sphere is opposite of that in the interior case (see the small ellipse 
outside the large sphere in Fig. [TTJa)). The continuation of the functions / and g 
from "interior" to "exterior" is displayed in Fig. [TT|a), where the transition from 
one orientation to the other is clear. 

Second, we compute numerically from Eq. ( [82) the interaction energy of a finite- 
size metal sphere with the cavity walls when the separation, d, between their sur- 
faces tends to zero. In this limit the Casimir force F between two conducting 
spheres, which is attractive, is proportional in magnitude to c/~ 3 , where d =R — r — a 
is the separation of surfaces. The coefficient of c/~ 3 is given by the PFA, 

r> n he rR 

\imd i F = -—— . (84) 

d^o 360 r + R 

This result holds for both the interior and the exterior configuration of two spheres. 
For fixed r we formally distinguish the cases: R > for the exterior, R — >• °° f or the 



plate-sphere, and R < for the interior configuration, see Fig. 11 b). All possible 
configurations are taken into account by considering — 1 < r/R < 1 . Although we 
know of no derivation of the functional form of the Casimir force beyond the leading 
term in the PFA, our numerical results are well fit by a power series in d/r, 
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7t 3 hc rR ( , , , d , , , d 2 \ 

We have used this functional form to extract the coefficient 6\ (r/R). 

Although the PFA is accurate only in the limit d/r — > 0, it can be extended in 
various ways to the whole range of d, r, and R. Depending on the surface O from 
which the normal distance to the other surface is measured, one obtains the "O- 
based" PFA energy. Clearly, the result depends on which object one chooses as O, 
but the various results do agree to leading order in d/r. We can choose either of 
the two spheres to arrive at the "r-based PFA" or the "R -based PFA", see Fig.fTTTb). 
Either one yields a 'correction' to the leading order PFA, 



where x = r/R. In Fig. 11 b) we plot the values of 9\ extracted from a numerical 
evaluation of the force from Eq. ( [82] i for various values of r/R < 0. For reference, 
the two PFA estimates are also shown. 

The numerical data in Fig. [TTTb) show a smooth transition from the interior to 
the exterior configuration. Although the PFA estimates do not describe the data, the 
r-based PFA has a similar functional form and divergence as x — > — 1 . Therefore, we 



fit the data in Fig. 11 b) to a function, 9i(x) = — (k\x + kix/{\ + x) + k^) and find, 
fci = 1.05±0.14,& 2 = 1.08±0.08,£ 3 = 1.38 ±0.06. Notice, however, that the actual 
function 9\ (x) is not known analytically and that the fit represents a reasonable 
choice which may not be unique. Our results show that the correction to the PFA 
has a significant dependence on ratio of curvatures of the two surfaces. 
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